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ABSTRACT 

We explore a possible origin for the puzzling anti-correlation between the formation epoch of 
galactic dark-matter haloes and their environment density. This correlation has been revealed 
from cosmological iV-body simulations and is in conflict with the Extended Press-Schechter 
model of halo clustering. Using similar simulations, we first quantify the straightforward as- 
sociation of an early formation epoch with a reduced mass growth rate at late times. We then 
find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects 
dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the 
assembly of haloes due to the shear along the large-scale filaments that feed the massive halo 
to actual mass loss in haloes that pass through the massive halo. Using the restricted three- 
body problem, we show that haloes are prone to tidal mass loss within 1 .5 virial radii of a 
larger halo. Our results suggest that the dependence of formation epoch on environment den- 
sity is a secondary effect induced by the enhanced density of haloes in filaments near massive 
haloes where the tides are strong. Our measures of assembly rate are particularly correlated 
with the tidal field at high redshifts z ~ 1. 

Key words: cosmology: theory, dark matter, large-scale structure of Universe - galaxies: 
formation, haloes - methods: N-body simulations 



1 INTRODUCTION 

Dark-matter haloes provide the gravitational potential wells in 
which galaxies form, and as such, the understanding of their as- 
sembly process during the cosmological history is a key element in 
the theory of galaxy formation. Focusing on gravitational physics 
while avoiding the complications associated with baryonic physics, 
the study of the growth of dark-matter haloes is the straightfor- 
ward, more solid part of the theory of galaxy formation. The pro- 
cess of bottom-up assembly of dark haloes is described successfully 
by an approximate formalism, the Extended Press-Schechter (EPS) 
or excursion-set model, which combines the statistics of an initial 
Gaussian random fluctuation field with the gravitational instability 
theory of linear fluctu ation growth and with the non-linear spher- 
ical collapse model jPress & Schechterl 1 19741 : iBond et al.l 1 19911 : 
lLacev & Colelll993l) . In the simplest version of this model, the as- 
sembly history of a halo is predicted to be fu lly determined by its 
mass, independently of its environment (e.g. iBardeen et al.lll98rj| : 
lEfstathiou et al.lll988l : lMo & Whitell996t) . More sophisticated ver- 
sions of EPS do involve a certain environment depend ence, but they 
do no t recover the correlatio n seen in the simulations iCatel arTet al.l 
1 19981 lTaruva&Sutoll200Cl) . The mere mass dependence, moti- 
vated by the simplest EPS model, lies at the basis of the widely 
used scheme for galaxy clustering that utilises the Halo Occupa- 
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tion Distribution (HOD) (see lCoorav & Shethl2002l and references 
therein). 

However, numerical simulations of cosmological structure 
formation revealed that the halo assembly hist ory does depend 
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Neistein & Dekell 120081: lAngulo et all 120081 : iLi et all 120081) . For 
example, when the halo history is parameterised by a formation 
redshift 2f rm corresponding to the time by which the halo had 
assembled one half of its current mass, haloes that formed earlier 
are found to be more strongly clustered and, equivalently, reside in 
denser environments. This effect is valid for haloes that are much 
smaller than the current characteristic gravitationally collapsing 
mass-scale M* (~ 10 13 Mq at redshift zero) and reverses for 
haloes above M*. The origin of this and related deviations 
from the predictions of EPS, sometimes termed "assembly bias" 
dGao & White 2003), is not obvious — it poses a very interesting 
open theoretical question which deserves a simple understanding. 

The re have been several attempts to understand the assem- 
bly bias. ISandvik et all d2007l) noted that the assembly bias can 
be incorporated in an excursion-set scenario where the spherical 
collapse model is replaced by ellipsoidal collapse and the density 
threshold for halo collapse depends on the shape and size of the 
local environment. Despite recovering some environment depen- 
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dence, the effect found by t hese authors is too w eak to explain the 
TV-body results. However, ISandvik et in £2007) considered only 
the tidal field generat ed by the linear (Gau ssian) density pertur- 
bations, as derived bv lDoroshkevicr] dl97dh . TV-body simulations 
indicate that at late times and small scales, the distributions of tidal 
field eigenvalues d epart significantly from those for a Gaussian ran- 
dom field (see e.g. lHahn et al.|[2007f) . Environmental e ffects in the 
ellipsoidal collapse model have also been discussed by Desjacques 
( 2008) but the model investigated by this author could not repro- 
duce fliecorrect sign of the assembly bias at low masses. 

IWang et alj [2007) found that small old haloes tend to reside 
next to very massive haloes. Considering the relative velocity of the 
dark-matter flow surrounding the haloes, they observed that parti- 
cles surrounding the haloes with highest Zf oxm have large relative 
velocities. 

Most recently, iDalal et all d2008h showed that the assembly 
bias at high masses follows naturally from the statistics of peaks 
of Gaussian random fluctuations when accounting also for the peak 
curvature. At low masses, these authors find a systematic trend of 
halo age w i th th e velocit y dispers i on of the environment. Both, 
IWang et all J2007l) and IDalal et alj d2008h suggest that the accre- 
tion onto a halo is suppressed because environmental velocities ex- 
ceed the halo's virial velocity. A possible scenario leading to these 
high velocities is that the dark-matter flu id is "hotter" in regions 
of higher density (see also lMo et all2005h after shell-crossing dur- 
ing the formation of the host large-scale structure. This is difficult 
to reconcile with other findings that the dark-mat ter flow in dense 
filaments is rather " cold" (e.g. iKlypin et alj|2003l) . In particular, 
ISandvik et al. I d2007l) argue that the mass of the large-scale envi- 
ronment in which the low-mass haloes form is not large enough 
to accoun t for this explanat ion. Furthermore , the m easures em- 
ployed by IWang et al] d2007l) and IDalal etail d2008l) are not sen- 
sitive to distinguish between various scenarios: "hot" flows, haloes 
that m ove against the flow of their environment on "unorthodox" 
orbits jLudlow et alj |2008) and strongly sheared flows. 

In this article, we address the hypothesis that the early forma- 
tion times of haloes of a given mass much below the non-linear 
mass in dense environments are primarily driven by the tidal sup- 
pression of the halo growth rate in the vicinity of a neighbouring 
massive hal o. Small, typical (low -sigma) haloes tend to form in fila- 
ments (e.g. lPorciani et alj2002al) . and stream along them into more 
massive, rare (high-sigma) haloes. The tidal field in these regions 
is likely to generate an ordered shear flow and thus suppress the 
halo growth by accretion or mergers. This suppression is expected 
to be correlated with the proximity to the massive halo. Those small 
haloes that have already passed through the inner regions of bigger 
haloes are likely to have lost mass due to tidal s tripping dGill et all 
l2005l : lDiemand et al.ll2007l : lLudlow et alj|2008h and thus show the 
earliest formation times. 

In order to test this hypothesis, we quantify in several differ- 
ent ways the halo assembly history and its environment, and study 
the correlations between the adopted quantities. The traditional pa- 
rameterisation of assembly history via the formation epoch Zform 
is complemented with more direct measures of halo growth rate 
at different redshifts. The environment, commonly referred to by 
the halo number density or the amplitude of the two-point corre- 
lation function, is characterised instead by the tidal/shear field, via 
the eigenvalues of the deformation tensor. The correlations between 
assembly history and environment are shown graphically and mea- 
sured via correlation coefficients. The primary role of the tidal field 
over the density field is to be addressed by the relative strength of 
the corresponding correlations with the assembly history. 



The outline of this article is as follows. In Section 2, we sum- 
marise the specifics of our TV-body simulations and the dark-matter 
halo catalogues. In Section 3, we recover the dependence of halo 
formation times on environment density and clustering amplitude. 
In Section 4 we introduce further measures of halo assembly com- 
plementing formation redshift. In Section 5, we define the relevant 
characteristics of the sheared flow around haloes and its tidal ori- 
gin before we illustrate its relation to mass assembly. In Section 6, 
we demonstrate the role of tidal effects in the assembly bias, by 
quantifying the correlation between the tidal field and the assembly 
rate, and comparing it to the correlation with environment density. 
We discuss our analysis in Section 7 and provide our conclusions 
in Section 8. 



2 NUMERICAL SIMULATIONS 

2.1 Specifics of the TV-body simulations 

We use a series of three high-resolution cosmological N-body sim- 
ulations thatwere obtained with the tree-PM code GADGET-2 
dSpringed r2005). These simulations are used to follow the non- 
linear evolution of density perturbations in a flat ACDM cosmology 
with a matter density parameter f2 m = 0.25, baryonic contribu- 
tion fib = 0.045, and a present-day value of the Hubble constant 
Ho = 100/ikms -1 with h — 0.73. The initial power spectrum 
has a long-wave spectral index of n = 1, and is normalised so that 
the rms fluctuation within a sphere of 8 /i _1 Mpc linearly extrapo- 
lated to the present time is as = 0.9. Each of the three simulations 
follows 512 3 collisionless dark matter particles in periodic boxes of 
sizes 45, 90 and 180 /i _1 Mpc, re spectively. Initial co nditions were 
generated using the GRAFIC tool dBertschingerl200 ll) at 2 ~ 79, 65 
and 52 for the three boxes. The corresponding particle masses are 
4.7 x 10 7 , 3.8 x 10 s and 3.0 x 10 9 ft _1 M Q , respectively. Particle 
positions and velocities were saved for 30 time-steps, logarithmi- 
cally spaced in expansion parameter a between 2 = 10 and z = 0. 
Unless otherwise stated, data from the simulation with the highest 
mass resolution, the 45 /t -1 Mpc box, is used for the analysis of the 
lowest mass haloes. Results have been checked for resolution and 
box size effects against the other boxes. 

2.2 Identification of haloes and the properties of their 
assembly 

We identify boun d structures using the standard friends-of-friends 
(FOF) algorithm jDavis et alj[l985h with a linking length equal to 
0.2 times the mean inter-particle distance. Among all these dark- 
matter haloes, we keep only those haloes at z = that fu lfil the 
virialization constraint \2T/V + 1 < 0.5 l lBett et al.ll2007l) where 
T is the kinetic energy and V is the potential energy of the group 
of particles. In addition, we also exclude those groups where the 
most-bound particle is further than J? max /4 away from the centre 
of mass, where 7? max is the distance of the furthest FOF particle 
from the group's centre of mass. 

For each halo at redshift z = we find its progenitors at any 
desired redshift 2 > by identifying all the haloes that have par- 
ticles in common with the final halo. At each snapshot, the most 
massive progenitor that contributes at least half of its particles to 
the final halo at 2 = is identified as the "main progenitor". We 
define the formation redshift 2f orm of a halo of mass Mo at 2 = 
as the redshift at which a main progenitor with mass M > Mo/2 
first appears in the simulation. The "assembly history" of a halo, 
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for our current purpose, is simply the mass-growth history of its 
main progenitor, M(z). 

In addition, we say that a halo underwent a "mass-loss" event 
if between any successive snapshots along its history it contributed 
at least 3% of its mass to a more massive halo that is not a main 
progenitor of the same final haloQ 



3 THE DEPENDENCE OF FORMATION TIME ON 
ENVIRONMENT 

The assembly bias has been detected using measures of either envi- 
ronment density or halo clustering. We first confirm these findings 
using our simulations. 

3.1 Dependence on Environmental Density 

At a fixed halo mass, the assembly bias can be measured as a sys- 
tematic trend of h alo formation time with the local density of the 
environment (e.g. lHarker et al]|2006h . We estimate the density at 
halo positions by smoothing the distribution of N-body particles on 
a given scale by convolving the particle positions with a spherical 
top-hat kernel of radius R. In practice, we first compute the density 
on a Cartesian grid using cloud-in-cell interpolation and then per- 
form the convolution with the kernel using FFT. We denote with the 
symbol Sr the mass overdensity computed using a smoothing ra- 
dius of Rh~ Mpc. Similarly, we call 8r 1 _ r 2 the overdensity com- 
puted within a spherical shell with an inner radius of R\ h~ 1 Mpc 
and an outer radius of R2 /i _1 Mpc. 

Figure [TJ (left panel) shows the median and the 16th and 84th 
percentiles of Zf as a function of S2-5 for ou r simulated haloes i n 
two mass bins. We recover the trend found by lHarker et alj d2006h . 
and extend it to lower masses. We note that the removal of mass- 
loosing haloes from the sample makes a difference preferentially 
at high densities and at the high end of the Zf distribution. We will 
return to these haloes in Section l631 where we quantify their abun- 
dance. 

3.2 Dependence on Clustering Amplitude 

Alternatively, the assembly bias can be discussed in terms of the 
halo clustering amplitude. The 2-point correlation function £ (r) is 
defined as the excess probability (with respect to random) to find 
two haloes in volume elements dVi and dV~2 separated by a co- 
moving distance r, 

dP 12 (r) = n 2 [1 + dVi dV 2 , (1) 

where n indicates the mean halo number density. The assem- 
bly bias has been originally seen as a systemati c dependence of 
the amplitude of £( r) with halo formation time dGao et al]|2005l ; 
iGao & Whitell2007l) . At constant halo mass, older haloes tend to 
be more strongly clustered than younger ones. The age dependence 
of the clustering amplitude is evident only for haloes with mass 
M < M*(z). Here M*(z) denotes the mass scale for which la 
density fluctuations typically collapse at redshift z according to the 
Press-Schechter model. The linear rms overdensity on this mass 
scale has then to be a(M,(z), z) = S c ~ 1.686. At all redshifts, 
the number density of haloes with M > M» (z) is exponentially 

1 This is equivalent to a ~ 5 per cent change in measured formation red- 
shift (see eq.[4}- 



suppressed. In Figure [TJ (right panel) we show how £(r) changes 
with halo age in our simulations. The ratio of the d ifferent corre- 
lation functions agrees very well with the results of IGao & White 
(120070 and extends the analysis also to lower masses. 



4 HALO MASS ASSEMBLY HISTORIES 

In this section we introduce measures of halo assembly com- 
plementing the common parameterisation of assembly histories 
through the formation time. 

4.1 The mass assembly rate of haloes 

We define the current assembly rate of each halo by taking the time 
derivative of its smoothed assembly history S[M(z)], where 5 is 
some smoothing operation. From the smoot hed assembly his tories, 
we define the current assembly rate as (cf. lNeistein et alj|2006l) 

A(z) = ~laS[M(K)]. (2) 

We apply smoothing by local regression (LOESS. Icievelandl 19791) . 
i.e. the data are fitted by a linear model inside a window and the fit 
value is evaluated at the centre of each window position. A window 
of 0.5 in redshift proved to be well suited and has been used in our 
analysis. It is necessary to smooth since the assembly histories of 
single haloes are prone to fluctuations that are due to the intrinsic 
stochasticity of the simulation/halo-identification process. 

In close analogy to the definition of the formation redshift, we 
will use the ratio G(z) = M(0)/M(z) between the final mass 
and the mass at redshift z as the "forward" measure of halo growth 
since redshift z. 

Both quantities A(z) and G(z) are strongly related with the 
integral measure «f orm (note that G(zform) = 2 by definition of 
Zform)- In Table[TJwe show the Spearman rank correlations between 
these different measures as a function of redshift for the main pro- 
genitors of haloes with masses between 2 and 4 x 10 10 and 2 and 
4 x 10 11 /i _1 Mq at z = 0. The median formation redshift for the 
first sample of haloes is 1.71 while it is 1.44 for the latter. We use 
the rank correlation coefficient rather than the linear Pearson cor- 
relation coefficient as it does not assume a linear relation between 
the correlated quantities. 

For haloes of a given mass at z — 0, a large (small) G(z) at 
high redshift implies trivially that the main progenitor mass was 
already high (low) at that epoch so that its formation time is high 
(low). The progenitor then also has a high (low) assembly rate A(z) 
at early times and a low (high) A(z) at late times. 

Thus, a high formation redshift is related to a low assembly 
rate A(z) at late times (2 < 2f or m) and a high assembly rate at 
early times (z > Zi OIla ). The latter can be seen for the higher mass 
sample where the median formation redshift is sufficiently low so 
that the correlation coefficient changes sign at z — 2. At late times, 
the formation redshift depends most sensitively on the assembly 
rate around 0.5 < z < 1 for the low mass sample and around 
2 ~ 0.5 for the higher mass sample. 

4.2 The influence of mass loss on formation time 

While typical analytic models of structure formation predict a 
monotonous growth of halo mass over time, haloes do undergo 
interactions with other haloes which can lead to tidal mass-loss. 
The influence of such a mass-loss event on the measured formation 
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Figure 1. (Left Panel): Median formation redshift (solid line) as a function of the overdensity £2-5 (see text for explanation) for all haloes, including those that 
experienced tidal mass loss (red), and only those haloes that never underwent such a mass loss (blue). The dashed lines indicate the 16th and 84th percentile 
in each bin. Results are shown for two different mass-ranges, indicated in the respective panels. The dotted lines show the median formation redshift for each 
mass bin. (Right Panel): Two-point correlation functions for all haloes (solid lines), and for those haloes of given mass that have 20 per cent highest/lowest 
formation times (filled squares/open circles). Errorbars correspond to Poisson errors. All results are shown including those haloes that experience tidal stripping 
of their mass (red) and excluding them (blue). See Section l2~2l for the definition of tidally stripped haloes. 
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Table 1. Spearman rank correlation coefficients p B between three different 
measures of halo mass assembly: the formation redshift Zf orm , tne assem- 
bly rate A(z) and the mass growth ratio G(z). The coefficients are given 
for the progenitors of haloes identified at 2 = with masses between 2 
and 4 X 10 10 h~ 1 M Q (top) in the 45 /i _1 Mpc box, and for haloes with 
masses between 2 and 4 X 10 11 /i _1 M (bottom) in the 90 h _1 Mpc box. 



time of a halo can be estimated as follows. S everal studies of struc- 
ture formation in a b ottom-up scenario (e.g. Ivan den Boschll200i 
IWechsler et alT2 002) have suggested that the bulk of haloes assem- 
bles the mass exponentially with redshift: 

M(z) = M(0) exp(-az) (3) 

where a is the mass-assembly rate. In this case, Zf olrn — a" 1 log 2. 
Let us assume that the smooth mass accretion history is interrupted 
at redshift z\ OBB by a single mass loss of A/i oaa = x M(z\ OSB ). For 
simplicity, the growth rate a is assumed to be identical before and 
after zi oss (it is trivial to extend our results to the most general case). 
A short calculation shows that mass loss leads to an increase in 



measured formation redshift given by 

1 2 

.zform = — In . (4) 

a 1 — x 



5 THE EFFECT OF TIDAL FORCES ON THE 
ASSEMBLY OF HALO MASS 

In this section we will demonstrate that the proximity to a massive 
halo provides a special environment for the formation and evolution 
of lower-mass haloes due to tidally sheared flows influencing their 
assembly. 

5.1 Halo evolution in the proximity to a massive halo 

For a large number of haloes in our simulation, we investigated 
the evolution of their Lagrangian patches and of their Lagrangian 
environments. At all epochs, the Lagrangian patch of a halo is de- 
fined by the set of particles that are identified as part of the halo 
at z = 0. Similarly, we define the Lagrangian environment as the 
set of particles that in the initial conditions of the simulations were 
lying within a cube of 3/i -1 Mpc centred on the corresponding La- 
grangian patch. We observe significant differences between old and 
young haloes of a given mass. 

The older haloes are preferentially found in pronounced large- 
scale structure already at high redshifts, while this is not true for the 
younger haloes. In Figuref2]we show the evolution of the particles 
making up two different haloes at z = together with their La- 
grangian environment. One of the haloes has a very high formation 
redshift of zf OIYn — 3.1 while the other halo has a low formation 
redshift of Zf orm = 0.42. It is evident from the evolution of the 
environment of the old halo that, once the host filament formed, 
its neighbouring haloes show an expanding, "Hubble-like", motion 
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Figure 2. Evolution of the Lagrangian patch and environment of two haloes with mass 1.4 X 10 11 ft _1 Mg at z = and formation redshift 2f orm = 3.1 
(left) and 2f orm = 0.42 (right). Particles that will constitute the halo at z = (the Lagrangian patch) are shown in red while particles that are contained in a 
3 7i Mpc cube centred on the patch (the Lagrangian environment) are shown in black. The centre of the coordinate system moves with the centre of mass of 
the particles that constitute the halo at z = 0. 



along the filament, which is not present for the younger halo. In the 
rest-frame of the halo, this expansion along the filament appears as 
a recession of the neighbouring haloes in that direction. Remark- 
ably, in both cases, the last major merger occurs perpendicular to 
the host large-scale structure. Thus, the gravitational shear acceler- 
ates the collapse perpendicular to th e filament while it slo ws down 
the collapse along the filament (cf. iPorciani et al.ll200 2a b). Tidal 
forces are expected to increase substantially in the vicinity of very 
massive haloes leading to a strong correlation between the tidal 
field and the density field on larger scales. This leads to the correla- 
tion of halo properties that are influenced by the external tidal field 
with the density field. 



5.2 Quantifying sheared flows 

The dynamics of dark-matter flows (and thus the accretion onto 
dark-matter haloes) is regulated by the gradient of the gravitational 
potential V<^>, its Laplacian V 2 <^> and the tidal tensor 



d 2 

didi 



(5) 



where 5y denotes the Kronecker symbol. The presence of a mas- 
sive halo that locally dominates the gravitational potential gener- 
ates a strong tidal field in its neighbourhood. Over time, the tidal 
field produces a tidally induced shear in the velocity field of the 
surrounding particles. The flow is stretched in the direction of 
the massive halo and compressed in the plane perpendicular to it. 



Since haloes pre ferentially reside in filaments or sheets (see e.g. 
Hahn et al. 2007), tidal forces tend to accumulate coherently. This 
should impact the mass accretion history of smaller haloes neigh- 
bouring the massive one. To probe the influence of tidal forces, we 
quantify the velocity field surrounding small-mass haloes by ex- 
panding the local (peculiar) flow with respect to the bulk motion of 
the halo up to linear order in the distance x from the halo centre 



dvi 



(6) 



We next perform a least squares fit of the nine components of the 
velocity gradient tensor within one and four virial radii around each 
halo (cf. Porciani et alj2002al lbT). It is then possible to write the ve- 
locity gradient tensor as the sum of a symmetric and an antisym- 
metric part 



dxj 



1 (<hk 

2 \dxj 



dvj 
dxi 



1 

+ 2 



dvi 
dxj 



dv 3 
dxi 



(7) 
(8) 



where Dij corresponds to the rate-of-strain tensor, and Qij de- 
scribes rigid rotation due to nonzero vorticity u>i = (V x v)j = 
—Cijk^ijk, where eyit is the Levi-Civita symbol. The rate-of-strain 
tensor can be further decomposed as Dij = 9 dij + Eij where 
6 denotes the velocity divergence (or expansion scalar) and Ejj 
is the traceless shear tensor with eigenvalues /ii ^ fj.2 ^ M3- 
For growing-mode fluctuations in the linear regime, Tij oc Ey 
and therefore the eigenvectors of these matrices are aligned. How- 
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0.5 



(\ 3 ,M/d 3 ) 0.62 
(At 3 ,M/d s ) 0.75 



0.64 
0.64 



0.61 
0.57 



Table 2. Spearman rank correlation coefficients p s between the largest tidal 
field eigenvalue due to the tidally dominant neighbour and the strain eigen- 
values A3 and shear eigenvalues ^3 , respectively. The coefficients are given 
for haloes identified at 2 = with masses between 2 and 4 X 10 10 li _1 Mp). 



ever, the proportionality of the shear and tide is broken by higher- 
order terms. In general, the eigenvectors of the tensors remain 
parallel but the eigenval ues evolve differently with time (e.g. 
iBertschinger & Jainl ll994i 2 l In fact, while tides act as sources for 
the velocity shear, the latter is also influenced by the velocity gra- 
dient, vorticity and the shear itself. In what follows, we use the 
velocity shear as a measure of the integrated effect of the gravita- 
tional tides after investigating the connection between the two in 
Sectionl53l 

Both random motion and nonlinear ordered motion will be su- 
perimposed onto the ordered linear deformation, as defined in eq. 
l(6j, so that we can define a velocity dispersion along the three main 
axes of the rate-of-strain tensor quantifying this residual by 

ah(((v-A,x- llJ xx).w,) 2 ), (9) 

where Ai A2 ^ A3 are the eigenvalues and Wj are the cor- 
responding eigenvectors of the rate-of-strain tensor Dij. Further- 
more, we define a simple relative velocity v rc i between a halo's 
centre of mass velocity and the mean velocity of the region between 
one and four virial radii. If the velocity field is completely linear 
then both |v re i| = and |er| = 0. We observe that, apart from 
very little scatter, \a\ is completely dominated by |v re i|. Nonzero 
|tr| and |v re i| can arise when the velocity grows non-linearly with 
distance from the halo, due to the presence of other haloes within 
1—4 virial radii or in special circumstances when the halo has in- 
deed a relative velocity with respect to its surrounding flow. We will 
return to this last aspect in Section [6731 Furthermore, noise in the 
velocity field due to finite numerical resolution will also increase 
\a\ by causing spurious deviations from a flow with constant strain 
and thus impacting the goodness of the fit (see also Section lrjTfl . 

5.3 The tidal origin of sheared flows 

In this section we will provide evidence that the tidal field due to 
a neighbouring massive halo induces shear in the dark matter flow 
surrounding smaller haloes in its vicinity. To this end, we will com- 
pare the constant strain approximation (CSA in what follows) in- 
troduced in the previous section, i.e. that the flow is described by 
Dij in eq. (8J, to a restricted three-body approximation for the tidal 
field. This will allow us to use the shear measured from the flow to 
quantify the importance of tides on the mass assembly of haloes in 
the rest of the paper. 

Consider the restricted three-body problem, consisting of a 
larger body of mass M , a smaller body of mass m C Af and a 
third of negligible mass. The third body can have stable circular 
orbits around the smaller mass m only within the Hill radius 



rmu 



Am) 



1/3 



(10) 



2 This holds in the Lag rangian frame before shell-crossing (see also 
Barnes & Row lingson 19^). 



where d is the distance b etween m and M (for a d etailed deriva- 
tion we refer the reader to lMurrav & Dermottll200ol) . Since in cos- 
mological context we do not have the idealized situation of the re- 
stricted three-body problem, we define the Hill radius of each halo 
as the minimal Hill radius due to any of its neighbour^]. Now, given 
that the eigenvalues of the tidal tensor Tij for a point mass M at 
distance d are given by 



,_ , _ 1 _GM_ _GM_ 2GM ) 

i T 'i = i d 3 ' d 3 ' —tf 



(ID 



this choice is equivalent to selecting that neighbour which has max- 
imal tidal influence GM/d z - the tidally dominant neighbour halo. 
The eigenvector t3 associated to the third eigenvalue T3 points ra- 
dially towards the mass M, while ti and t2 span the perpendicular 
plane but are otherwise undetermined due to the symmetry of the 
problem. 

By investigating the flux through the Hill sphere (with volume 
VHiii and surface 9VHiii), it is then possible to connect the sheared 
motion of the flow around a halo, measured using the CSA defined 
in Section l572l to the tidal forces due to the neighbouring halo. This 
flux is given by 



4> 



VHiii 



vdA 



47T 

V ■ v dV ~ y-mn (Ai + A 2 + A 3 ) 



(12) 



(13) 



(A1 + A2 + A3), 



(14) 



where the first equality is just an application of the divergence the- 
orem and the second follows after applying the CSA to the flow 
where the A^ are assumed to be constant within rnm. 

In order to verify that the shear measured in the flow indeed 
satisfies the CSA, we can rewrite eq. d!3t as 

(Vr) „ 1 
J-ffill 3 

where (v r ) is the mean radial velocity at rmu. Written this way, 
the left hand side is independent of the assumption of constant 
strain. We can thus verify the validity of equation d 14b directly in 
the A^-body simulations in order to assert that the CSA is justi- 
fied. In Figure [3] we show (v r ) /rnm as a function of the mean 
of the strain eigenvalues for haloes with rum < 4r v i r , the range 
used when fitting the CSA flow model. The mean radial veloc- 
ity was measured using the closest 100 particles outside the Hill 
sphere. Results are shown for haloes that have final masses at 
z — between 2 — 4 x 10 /! _1 M . The rank correlation is 
Ps({v r ) / r mii,X^ Ai/3) = 0.74. This very strong correlation in- 
dicates that the constant strain model serves well to predict the flux 
into the Hill sphere. 

3 This will possibly overestimate the Hill radius. The restricted three-body 
problem applies to the orbit of a test particle of zero mass around the mass 
m. While in the case of smooth accretion the third body can indeed be 
thought of as having negligible mass, in the case of mergers, the restricted 
three-body problem is of course of limited applicability. 

Furthermore, the presence of finite masses other than m and M makes it 
unfortunately impossible to use radial acceleration vectors to find the region 
of gravitational influence of the second body m. Any self-bound mass (e.g. 
subhaloes or other tiny haloes outside the virial radius) will produce a dip 
in the potential. In order to use acceleration vectors it would be necessary 
to remove all structure smaller than the halo itself, which is practically im- 
possible. In fact, for the same reasons we consider the shear tensor rather 
than the tidal tensor since the latter is affected by the same problems. Shear 
is the integrated result of tidal forces and thus probes directly their effect on 
the flow at the relevant scales. 
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Figure 3. Comparisons between the flux onto haloes from particle velocities and in the CSA model, and between the tides due to the tidally dominant neighbour 
halo and measured shear in the CSA. (Left panel:) The flow through the Hill sphere A;/3 from the linear fit to the flow against median (v r ) /rjjm measured 
directly from the simulation. Results are shown for haloes with r Hi n ^ 4r v ; r , the region used for determining the A^. (Right panel:) Median eigenvalue of the 
shear tensor [13 measured from the flow as a function of the tidal force due to the tidally dominant neighbour halo T3 = 2GM/R 3 . Errorbars show the 16th 
and 84th percentiles in each bin. Bins are chosen to include equal numbers of datapoints. The dotted lines indicate x = y. Data are given for the progenitors 
of haloes with masses between 2 and 4 X 10 10 h^M®. 
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Figure 4. The probability distributions of angles | cos t9 4 1 between the three 
rate-of-strain eigenvectors Wj and the direction to the tidally dominant 
neighbour halo for the progenitors of haloes with 2 = masses between 2 
and 4 X 10 10 /i~ 1 M for three different redshifts: 2 = (blue), 2 = 0.5 
(green) and 2 = 1 (red). The black dotted line indicates a random distribu- 
tion of angles. 



In addition, we can check directly whether A3 is probing the 
tidal influence due to a neighbouring halo in which case we expect 
that A3 oc T3. Since T3 oc M/d 3 is a positive definite quantity, we 



show the shear eigenvalue ^3, which has the same property, as a 
function of M/d 3 in the right panel of Figure [5] Results are given 
at three redshifts for the progenitors of haloes with masses between 
2 and 4 x 10 10 ft. _1 M . We find a tight connection between the 
two quantities which is also reflected in the high rank correlations 
given in Table [2] 

Finally, we can also verify whether the tidal field eigenvector 
t3, which is simply the normalized vector pointing radially towards 
the tidally dominant halo, is parallel to the third eigenvector W3 of 
the strain tensor and perpendicular to the first and the second, as 
would be expected if tides are inducing shear in the flow. Figure 
[4] shows the distribution of the cosine of the angles between these 
vectors, i.e. 



I cos OA 



Wi ■ t 3 



(15) 



at three redshifts in the simulation. If the strain on the fluid is in- 
deed caused by the tidal field due to the neighbouring halo, then 
an alignment angle | cos Oz \ — lis expected which is consistent 
with the results shown in Figure [4] It is interesting to see that the 
degree of alignment - reflected in the width of the distributions - is 
slightly decreasing with decreasing redshift. 

The findings presented in this section - in particular (1) that 
the rate-of-strain eigenvalue A3 is very strongly correlated with the 
tidal field eigenvalue due to the tidally dominant neighbour halo 
and (2) that the associated rate-of-strain eigenvector points into the 
direction of this halo - lead us to conclude that A3 can be used 
as a proxy to quantify the integrated effect of tides on the flow 
surrounding haloes in our subsequent analysis. 
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Figure 5. (Left) Cumulative distribution functions (CDF) of Hill radii in units of virial radii, and (Middle) CDFs of the mass flux through the Hill sphere, 
(Right) same as middle panel but only for low-velocity particles with v < t> v ir- All data are shown for the progenitors at z = 0.5 of haloes with final masses 
between 2 and 4 X 10 10 h~ 1 Mq (and also at z = for the left panel). The CDFs are shown in black for the sample of all these haloes and in red (blue) for 
those haloes with the 20 per cent lowest (highest) assembly rate A(z). Stripped haloes are excluded. 



5.4 The tidal influence on mass assembly 

The assembly of haloes is regulated by the amount of accretable 
material in their vicinity. The sphere of gravitational influence of 
a halo is bounded by the Hill sphere, only particles within it are 
potentially accretable. We will in this section quantify differences 
in the supply of accretable material and its relation to mass growth 
and tides. 

In Figure|5] left panel, we show the cumulative distribution of 
Hill radii (in units of virial radii) for the z = 0.5 progenitors of 
haloes with masses between 2 and 4 x 10 10 /i -1 Mq and for these 
haloes at z=0. The data are shown for all haloes and for those with 
the 20 per cent highest and lowest current assembly rates A(z). We 
find that at both redshifts the estimated Hill radii of the least accret- 
ing haloes are significantly smaller than those of the most accreting 
haloes, which are comparable to those of all haloes. This implies 
that haloes with the lowest assembly rates experience stronger tidal 
fields than all other haloes. 

The mass assembly rate A(z), which can also be thought of as 
the mass flux through the halo surface 

dMh f a* 

— — = pv dS, (16) 

is expected to be correlated through the continuity equation to 
the supply rate of accretable material. The latter is the mass flux 
through the Hill sphere 

4>m\l= [ pv-dS. (17) 

In both definitions p is the matter density, is the halo mass. This 
correlation is however likely to be weakened by several factors. A 
particle does not enter the Hill sphere and the virial radius simulta- 
neously, especially when the Hill radius is large, leading to a time 
delayed correlation with the delay depending on the size of rum. 
Furthermore, the Hill radius shrinks linearly with the distance to a 
more massive halo (see eg. 1101 ). In addition, the tidal field will per- 
turb particle trajectories also inside the Hill sphere. Thus, entering 
the Hill sphere is a necessary but not sufficient condition for a par- 
ticle to be accreted in the future. Nevertheless, we can compare the 
flux through the Hill sphere for haloes with high and low assembly 
rates to probe the supply of accretable material and its relation to 
the assembly behaviour of haloes. In Figure [5] middle panel, we 



show cumulative distributions of the flux through the Hill sphere 
for haloes with highest and lowest current assembly rates A(z). 
The flux was computed using the closest 100 particles outside the 
Hill sphere from 

0Hiii^47r7?| m (v r p), (18) 

where v r is the radial velocity and (■) denotes the average over the 
shell occupied by the 100 particles. We find that the Hill sphere 
flux of the least accreting haloes is shifted towards more positive 
values relative to the most accreting sample, i.e. the flux into the 
Hill sphere is reduced for the least accreting haloes. In particular, 
these haloes have a median flux very close to zero. 

Particles that enter the Hill sphere might not be accretable if 
they have too hig h velocities to be captured by the halo (see also 
IWang et alj2007l) . In the plane perpendicular to the eigenvector be- 
longing to T3 the tidal field compresses the flow and thus acceler- 
ates particles onto the halo. Furthermore, high velocities could also 
be due to particles orbiting the bigger halo which has a larger cir- 
cular velocity than the smaller halo. To quantify the abundance of 
these particles that are potentially too fast for the halo to be cap- 
tured, we determine the fraction of particles entering the Hill sphere 
that have a velocity smaller than the circular velocity. This is a con- 
servative requirement since the halo should be able to capture par- 
ticles up to its escape velocity (which is twice its circular velocity). 
We thus recompute the flux through the Hill sphere using only these 
accretable particles. In Figure|5] right panel, the cumulative distri- 
butions of accretable fluxes are shown for the same halo samples 
as before. We find that the difference between the distributions in- 
creases slightly, the median flux through the Hill sphere is still very 
close to zero (if not slightly positive) for the least accreting haloes 
while it is negative for the most accreting haloes. 

These results suggest that accretion is not primarily limited 
by high velocities since the difference between the flux distribu- 
tions does not change substantially when considering all particles 
as opposed to only slow particles. Hence, non-accreting haloes do 
not even draw on a supply of slow particles suggesting that they 
are no longer convergence points in the local flow of dark matter 
particlesQ. 

4 Note that the relatively low flux through the Hill sphere for the most- 
accreting haloes appears to be in conflict with their assembly rates and 
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Figure 6. Cumulative distribution functions of the eigenvalues A^ of the 
deformation tensor -D;j. The data is shown for the z = 0.5 progenitors of 
haloes with masses between 2 X 10 10 and 4 X 10 10 h _1 MQ at z = 0. 
Lines correspond to all these haloes (solid, black), those 20 per cent with 
highest current assembly rate (short dashes, blue) and those with the lowest 
(long dashes, red). Stripped haloes are excluded. 



5.5 The shear flow near early forming haloes 

We will now discuss the connection between a sheared velocity 
field surrounding low mass haloes and their mass assembly prop- 
erties, thus providing further evidence for an influence of external 
tides on mass assembly. 

In Figure [6] we show the cumulative distribution functions of 
the eigenvalues Ai of the deformation tensor Dij for haloes with 
highest and lowest assembly rates A at z — 0.5, where the effect is 



thus the growth they should experience. This incongruity is however eas- 
ily resolved by considering the relative sizes of the Hill radii (see Figure 
f5). The flux is constrained to the assembly rate at the virial radius and to 

fa ~ -32npH(z)r 3 v . r ~ -9 x 10 11 (y%) 3 A^Ms/Gyr at the 
turnaround radius (assuming that the density at the turnaround radius is p, 
that r ta = 2r v j r and velocity is —H(z)r ta ). For the haloes we consider, 
r vir ~ 0.09 7l _1 ]vlpc and so </> la ~ -7 x 10 s h- 1 M e /Gyr. Thus, if all 
the material that the halo will accrete until z = is already within its 
Hill sphere, then the Hill radius will be on the order of or larger than the 
turnaround radius. Given that the overdensity within the turnaround radius 
in the spherical collapse model is found to be 4.56, we can simply com- 
pare the overdensity within the estimated Hill radius with this number to 
compare their relative sizes. We find that the median overdensities within 
the Hill sphere for the haloes with the 20 per cent highest and lowest as- 
sembly rate A(z) at z = 0.5 are 5.54 and 12.72, respectively. Thus, for 
the high-A(z) sample, the relatively low Hill flux is perfectly compatible 
with their large assembly rates as this is the flux determined close to the 
turnaround radius. Since their estimated Hill radii are comparable to their 
turnaround radii, the growth of these haloes is simply unaffected by tides. 
Indeed, 94 per cent of the high-A(z) and 99 per cent of the low-yl(z) sam- 
ple have at least their final mass contained inside their estimated Hill radius 
at z = 0.5. Note also that in dark-energy cosmologies the value for the 
turnaround overdensity at low redshifts can be significantly higher than the 
fiducial Einstein-de Sitter value (see e.g. lHorellou & BergeH2005f) . 



Figure 7. Median velocity field along the three main axes of the rate-of- 
strain tensor wi (top), W2 (middle) and W3 (bottom) in units of the halo 
circular velocity v c . The data is shown for the z = 1 progenitors of haloes 
with masses between 2 X 10 10 and 4 X 10 10 /i _1 M Q at z = 0. The sample 
of haloes is split into those 20 per cent with largest (blue) and smallest (red) 
increase in mass. Stripped haloes are excluded. Errorbars correspond to the 
medians (over all haloes) of the 16th and 84th percentiles (computed for 
each halo). The vertical dashed lines indicate the median Hill radii for the 
two samples of haloes. 



strongest. We observe a systematic shift of A2 and A3 towards more 
positive values for those haloes that grow slower. This implies that 
the flow around these haloes is more strongly sheared. As discussed 
in Section |531 sheared motion in the flow provides clear evidence 
for the presence of an external tidal field. In particular, for the low- 
A haloes, over 80 per cent have A3 > 0, i.e. the velocity component 
along the eigenvector associated with A3 is positive, particles are 
moving away from the halo in that direction. 

In Figure [7j we show the median velocity of the environment 
as a function of distance from the halo centre projected along the 
three main axes of the rate-of-strain tensor (which coincide with the 
main axes of the shear tensor). Results are shown for haloes grow- 
ing most and least in mass between redshift one and zero for the 
progenitors of haloes of fixed mass at z = 0. The velocity is given 
in units of the circular velocity of each halo v c = \J GM/R v i r . 
While the flow fields are comparable along wi and W2 apart from 
slightly larger scatter for the less growing haloes, these haloes show 
a qualitatively different flow along the eigenvector W3 correspond- 
ing to the largest eigenvalue A3. For the least growing haloes, along 
this direction, matter is receding from the halo centre leading to a 
Hubble-like flow that counteracts infall onto the halo. The data are 
shown at z = 1 where the difference in the flow field is strongest 
when splitting the halo sample by growth G(z). The scatter shown 
in Figure [7] indicates that there is a non-negligible fraction of par- 
ticles surrounding the halo that have velocities on the order of the 
halo's circular velocity, as discussed and quantified in the previ- 
ous section. Note that outflow along one direction is equivalent to a 
non-convergent velocity field. Such a non-convergent flow is in ex- 
cellent agreement with other detectio ns of matter outflo ws around 
haloes below the non-linear mass (e.g lPrada et al]2006h . 
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Table 3. Spearman rank correlation coefficients p s between measures of 
halo mass assembly and measures of environment over redshift. Measures 
of environment are A3 , the largest eigenvalues of the rate-of-strain tensor, 
the overdensity 8$ in the shell within one and four virial radii around each 
halo, 82 the matter field smoothed with a spherical top-hat filter of scale 
2 /i -1 Mpc and the deviation from a linear flow \er\. Measures of mass as- 
sembly are the formation redshift Zf orm , the assembly rate A(z) and the 
mass ratio G(z). The coefficients are given for haloes with masses between 
2 and 4 X 10 10 h~ 1 Mq at 2 = 0. These haloes have a median forma- 
tion redshift 2 mc( j = 1.71. The median mass Af mc d(z) of the halo sample 
at redshift 2 is given in units of I0 10 h~ 'Mq. Data is obtained from the 
45 h~ x Mpc box. Tidally stripped haloes are excluded. These results illus- 
trate that: (1) mass assembly is modulated by A3 and (2) the correlation of 
A3 with density grows over time. 



6 THE ROLE OF TIDAL EFFECTS IN THE ASSEMBLY 
BIAS 

The hierarchical structure-formation picture based on the Extended 
Press-Schechter model predicts that the mass of the main progeni- 
tor of a halo is determined by the sequence of "density peaks" on 
increasingly larger scales that exceed the threshold for collapse. In 
this picture, formation time depends only on the halo mass without 
any influence of the environment. 

We will argue in this section that in TV-body simulations the as- 
sembly behaviour of haloes is modulated by tidal effects ranging 
from the suppression of halo growth to tidal mass loss due to an 
encounter with a massive halo in the most extreme cases. 

6.1 The influence of tides on halo mass assembly 

In order to statistically quantify the interplay between environment 
and halo mass assembly, we present Spearman rank correlation co- 
efficients between measures of environment and measures of halo 
mass assembly in Tables[3]and|4]at three redshifts for the main pro- 
genitors of haloes identified at 2 = 0. In Table[3] data is shown for 
haloes in the 45 ft" 1 Mpc simulation with masses between 2 and 
4 x 10 10 h~ 1 M Q at z = 0. Table|4]shows the corresponding data 
for haloes with masses between 2 and 4 x 10 11 h~ 1 M ( v ) at z — 
from the 90 ft -1 Mpc simulation. 

We find that for both mass ranges the correlation of formation 
redshift and the eigenvalue A3 is strongest at redshifts close to the 
median formation redshift of the sample. Furthermore, the corre- 
lation of assembly rate A(z) with A3 increases between redshifts 
1 and 0.5 and is significantly stronger for the lower mass sample. 
Finally, we find that at all redshifts a large A3 is strongly correlated 
with a small growth G(z) and vice versa. 
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Table 4. Same as Table fj] but for haloes with masses between 2 and 4 X 
10 11 /i -1 Mq at 2 = 0. These haloes have a median formation redshift 
2 mc( | = 1.44. Data is obtained from the 90 ?t~ 1 Mpc box. 
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Table 5. Z-scores Z(x, y, 2) reflecting the significance in units of a that 
p s (x,z) > p s (y, z) for the correlation coefficients given in Tables[3]and|4] 
Negative values indicate the score for the alternative hypothesis p s (x, 2) < 
Ps(y,z). 

The correlations with A3 can now be compared with the de- 
pendence of assembly on environment density. At all redshifts 
above zero, the correlation coefficients between the overdensity in 
spheres of 2 ft _1 Mpc and the three measures of mass assembly are 
significantly weaker than those with A3 . Most remarkably, the cor- 
relation between 82 and 2f orm is strongest at z ~ 1. However, at 
these times «f or m correlates significantly stronger with A3 than with 
82 indicating that the density itself is not the primary cause for the 
suppression. As expected, the correlation between density and A3 
grows strongly with time. 

The significance of the quoted correlations can be assessed by 
testing the hypothesis that p s (x, z) > p s (y, z) for three quantities 
x,y,z respecting the fact that also x and y are corre lated. This 
can be calculated using the method of lMeng et all h992tPI . The Z- 
scores (significance in units of a) for the correlations in Tables [3] 
and [4] are given in Table [5] At redshifts 0.5 and 1, for all quoted 
correlation coefficients between A, G and Zf OIln with A3 and 82, 
the differences are significant at more than 3<r. 

5 The original iMeng et all jl992h test applies for Pearson correlation co- 
efficients but can be easily adopted to rank correlations by replacing all cor- 
relation coefficients and multiplying the standard error of the z-transformed 
values w ith 1.03 (this is valid for samples larger than 10 and p s < 0.9, see 
IZaj2007l) . 
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Figure 8. Two-point correlation functions for haloes with masses between 2 
and 4 X 10 10 at redshift zero. The two panels correspond to these haloes at 
redshift zero (top) and their main progenitor haloes at redshift one (bottom). 
The samples have been split into those haloes at the respective redshift with 
the 20 percent highest/lowest values in A3, the largest eigenvalue of the 
rate-of-strain tensor. Tidally stripped haloes are not included. Haloes ex- 
periencing strong tides, as probed by A3 , are much stronger clustered than 
haloes in low shear regions. 

We will now briefly discuss the influence of numerical resolu- 
tion effects in the simulation on these results. Since the mass range 
of haloes shown in Table[4]is resolved with 8 times more particles 
in the 45 ft -1 Mpc box, albeit for a smaller sample, we compared 
all correlation coefficients in the two simulations and found no sig- 
nificant deviations except for the fit residuals <x|. At all redshifts, 
correlations of assembly with |o"| were weaker and those with A3 
slightly stronger in the simulation with higher mass resolution. The 
observation that the fit residuals depend on resolution is indicative 
of a minor resolution dependence in determining the tensor eigen- 
values from the flow that results in a slight weakening of correla- 
tions with shear. Note also that the quoted results do not change 
significantly if haloes are selected in a fixed mass range at some 
redshift z > instead of at z — except that all correlations with 
|tr| are significantly weaker at redshifts z > again reflecting a 
resolution dependence of the goodness of our CSA model fit to the 
flow. 

6.2 Clustering behaviour of haloes in sheared regions 

Tidal fields are strongest in the vicinity of massive haloes where 
also the clustering of smaller haloes is substantially enhanced. In 
Figure [8] we show the correlation functions for haloes with the 
largest/smallest eigenvalue A3 of the rate-of-strain tensor at red- 
shifts zero and their main progenitor haloes at redshift one. As ex- 
pected, haloes at redshift zero with high A3 are significantly more 
clustered than average haloes while those with low A3 are less 
clustered. At redshift one the situation is slightly different: only 
haloes with highest A3 are significantly more clustered than average 
haloes. Furthermore, haloes with the lowest A3 show a strong atten- 
uation of the correlation function at distances below ~ 1 fa -1 Mpc. 



Figure 9. Fraction of haloes that have rniii < r v i T at some point during 
their evolution (black symbols) and fraction of haloes loosing at least 3 per 
cent of their mass to a more massive halo (grey symbols) as a function of 
halo mass at redshift z=0. Errorbars indicate the la Poisson uncertainties. 
Results are shown for the three simulations of different resolution described 
in the text. Dotted lines show exponential fits to the data points. 

This is a clear indication that these haloes are only found in more 
underdense regions. See also the discussion in Appendix lA2l for the 
dependence of clustering amplitudes on environmental overdensity 
in the initial density field. 

The tidal field, as quantified by the sheared flow around a halo, 
thus shows a strong correlation with all quantities of mass assembly 
as well as the dependence on clustering amplitude that is necessary 
to account for the halo assembly bias. 

6.3 Tidal stripping of halo mass 

A halo is prone to tidal mass loss once its Hill radius run, defined 
in eq. jlOh with respect to a massive neighbour becomes smaller 
than its virial radius r v i r . Defining haloes as objects with a fixed 
overdensity A v j r , we have m oc A v i r r^ ir for the smaller halo and 
M oc A v i r i?^ ir for the larger halo. Then the condition that run < 
r v ir is equivalent to 

d < 3 1/3 i?vir ^ 1.5-Rvir. (19) 

Thus, a halo is likely to suffer tidal mass loss once it has approached 
a more massive halo to within 1.5 virial radii, assuming that the 
restricted three-body approximation is justified and that the small 
halo remains nearly spherical. Haloes that tidally loose mass to an- 
other halo should be considered as the most extreme case of tidally 
suppressed mass growth: the range of their gravitational influence 
has become smaller than their virial radius. 

We will now quantify the population of these extreme exam- 
ples of tidally truncated mass assembly. As described in Section 
14.21 the stripping implies an underestimation of halo mass at 2 = 
relative to the assembly history at higher z, which is equivalent to 
an overestimation of the formation redshift. It is thus important for 
our analysis to quantify the abundance of these objects. 

In Figure [9] we show in grey the fraction of well resolved 
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haloes in our simulations that are found as isolated haloes at z = 
and lost at least 3 per cent of their particles to a more massive halo 
(Note that this definition entails haloes that are actually detected as 
subhaloes at z > 0). At a given mass, this fraction is on the order of 
a few per cent. The fraction is higher for smaller halo masses and 
reaches a value of ~ 4.5 per cent for M ~ 2 X 10 10 ft _1 M Q . 

However, these haloes do not represent the entire population 
of haloes that underwent mass-loss in an encounter with a more 
massive halo. First, the time resolution of our saved simulation 
snapshots does not permit to detect all haloes that were at some 
point subhaloes of another halo, and second, it is not guaranteed 
that the tidally unbound mass will be part of the more massive halo 
after stripping. In order to get an upper limit on their abundance, we 
consider now all those haloes for which the Hill radius is smaller 
than the virial radius at some point during their evolution and that 
are thus prone to tidal mass-loss. The fraction of haloes fulfilling 
this criterion is on the order of 20 per cent at the lowest masses 
we considered (see Figure|9] black line an d sy mbols) in ag r eemen t 
with the results of IDiemand et al. I J2007h and iDalal etafl 12008). 
The lower number of those haloes for which actual mass-loss is de- 
tected implies that the tidally unbound mass typically does not end 
up in the bigger halo. 

In Figure Q] (left panel) we showed how the correlation of for- 
mation redshift with overdensity 62-5 changes when the haloes 
that lost mass due to an encounter with a more massive object 
are removed from the sample. The rank correlation decreases from 
ps (62-5, zform) ~ 0.18 for the entire sample to 0.12 after exclu- 
sion of the stripped subset. Removing also the larger sample of 
those haloes that had runt < r v i r during their evolution results 
in a slight increase of the correlation coefficient to 0.16 (all num- 
bers are slightly higher when using S2 instead of 82-5). Note that 
the criterion rum < »"vir is roughly equivalent to selecting those 
haloes that were within 1.5 virial radii of another halo. Our results 
thus suggest that tidally stripped haloes are only the most extreme 
cases of tidally truncated mass assembly. 

Furthermore, haloes that are at the apocenter of their orbit af- 
ter passing through a massive halo are also likely to have very large 
relative velocities with respect to their environment. Considering 
the relative velocity between a halo and the surrounding region 
within 1 and 4 r v ir (cf. Section [53t . we find that haloes identified as 
stripped haloes indeed show large relative velocities. These relative 
velocities will lead to extreme velocity d ispersions a r ound haloes 
and are possibly dominating the results of lWang et al.l |2007) (their 
Figure 8, where median relative velocities are much larger than the 
halo virial velocity only in the tail of the formation redshift distri- 
bution). When removing the stripped haloes, we find that the cor- 
relation between relative velocity and formation redshift weakens 
substantially. 



7 DISCUSSION 

The formation redshift of dark-matter haloes is highest for haloes 
that have the lowest assembly rate at late times (cf. Section |4~TV Ty- 
ing this suppression of late growth to a physical effect that is most 
effective in dense/strongly clustered regions can thus potentially 
explain the correlation of halo clustering and formation redshift, 
the so-called "assembly bias". 

Tidal forces induce a sheared flow in the vicinity of small 
haloes that flow along filaments of the cosmic web (Section |5T| l. To 
quantify this effect, we measured the local influence of the larger 
scale tidal field by considering the rate-of-strain tensor of the flow 



in the immediate vicinity of haloes. We found that the progeni- 
tors of haloes of a given mass at z — that have the lowest mass 
assembly rates experience significantly stronger shear during their 
assembly history than those haloes with high mass assembly rates. 

Tidal effects are expected to be strongest close to massive 
haloes. Whenever a mass M at distance r from a smaller halo 
completely dominates the potential, the tide in the direction of this 
mass is given by 2GM/r s , while the perpendicular tides are given 
by —GM/r 3 . In the rest-frame of the smaller halo, matter is re- 
ceding along the direction to the massive halo. Furthermore, upon 
approaching the massive halo, the environment density will grow 
substantially and as a consequence the correlation between density 
and shear grows with time (see Section |6j. Since the abundance 
of small haloes is enhanced near massive haloes, the population of 
tidally suppressed haloes will be biased such that a dependence of 
formation redshift on clustering amplitude and environment den- 
sity arises as a secondary effect. The more fundamental character 
of tides is reflected in the stronger correlation of measures of halo 
growth with A3 than with environment density (Section|6j. 

The tidal influence of a larger halo on the accretion proper- 
ties of neighbouring smaller haloes can be readily understood from 
the restricted three-body-problem. Stable orbits around the smaller 
halo exist only within its Hill sphere (see Section [53V A halo can 
only accrete particles (and other smaller haloes) in the presence of 
a massive halo if they are contained in the Hill sphere. The radius 
of the Hill sphere shrinks linearly with the distance to the massive 
neighbour so that the amount of matter available for accretion is 
becoming increasingly smaller. In addition, a halo is prone to mass 
loss when the Hill radius becomes smaller than its virial radius. 
As shown in Section l63l this happens roughly at a distance of 1.5 
virial radii of the larger halo. 

In the most extreme cases, small haloes can undergo very 
strong tidal mass-loss when passing through the massive halo. A 
minor fraction of haloes of a given mass are indeed not typical 
isolated haloes since they are once again found outside the virial 
radius of a massive halo after their first pericenter passage through 
its potential well. These haloes ha ve also been observed in resim- 
ulations of single massiv e haloes dGill et alll2005l ; IDiemand et al.l 
l2007l ; lLudlow et alj|2008h . 

While tidal mass-loss inside the virial radius of another halo 
is the most violent and non-linear example, tidal suppression of 
halo growth seems to be alr eady present in the quasi-linear regime. 
iKeselman & Nus ser (2007) recover the assembly bias also using 
the punctuated Zel'dovich approximation where the mass of haloes 
grows monotonously by definition and tidal mass-loss is thus im- 
possible. 

We explore the link between assembly bias and the initial lin- 
ear density field in Appendix lAl where we compare the masses of 
haloes measured from the simulations with the masses they should 
reach in an isolated environment according to the EPS model. Our 
results indicate that the EPS predicted masses are roughly a factor 
of 2-3 higher for the oldest haloes than for the youngest. Further- 
more, we find a significant correlation between this mismatch and 
both the strength of the shear and large scale density at z — in 
the simulation. This provides further evidence that indeed the en- 
vironment, that is not considered when determining the EPS-like 
masses, prevents the mass to grow to the predicted value. Hence, 
while EPS would predict that the halo will assemble mass, the 
strong shear that EPS does not account for prevents this (see also 
van de We ygaert & Babulll994l) . If the correct final mass is known, 
it can be shown that the density dependence of formation redshifts 
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can be reproduced using a simple peak-background-split model 
(see Appendix [At. 

A question that is opened by our analysis is whether the last 
major merger of small haloes occurs preferentially perpendicular to 
their host filament such that the mechanism we described is poten- 
tially also related to alignments between ha lo spins and shapes and 
the surrounding large-scale structure (e.g |Pekellll985h . For small 
galactic dark matter haloes such strong alignments between halo 
shape and the direc tion to the nearby cluster have been found by 
iPereira et alj j2008h extending out to many virial radii of the clus- 
ter and by alignments of both ha l o shap e an d spin with host fil- 
aments and sheets by Hahn et al. (2007) and Aragon-Calvo et al. 
(2007) (see also[B asilakos et alJl2006l : lRagone-Figueroa & Plionisl 
l2007t) .[L ee & Er dogdu ( 2007) claim observational evidence for an 



alignment of galaxy spins with the reconstructed tidal field. In the 
picture outlined in this paper, the decreased probability for merg- 
ers or accretion along the filament would preserve the spin vector 
alignment with the filament/sheet which is established when merg- 
ing occurs preferentially perpendicular to the filament/sheet. 

We focused in this paper on the role of tidal effects in 
the assembly of haloes in different environments. We note that 
other global properties of haloes, such as the substructure in 
them and their spin parameter, also corre l ate with environ ment 
dWechsler et ail |2006| - iGao & White! 120071; iBett et alj l200l see 
e.g.). The possible relation between these properties and the sup- 
pression of accretion onto haloes should be explored. 



8 SUMMARY AND CONCLUSIONS 

Using a series of cosmological iV-body simulations, we demon- 
strated that the mass assembly of haloes is influenced by tidal ef- 
fects in the vicinity of massive haloes. The resulting suppression of 
mass growth at late times is reflected in an increase in formation 
redshift. The space density of small haloes near massive haloes is 
enhanced (since their spatial cross-correlation is high and positive). 
Thus, as a secondary effect, a tidally suppressed assembly of halo 
mass directly translates into biased spatial distributions of young 
and old haloes identified at z = 0: old haloes will be found more 
likely in regions of high density than younger haloes. 

The presence of tides has an important effect on the infall be- 
haviour onto a halo. To quantify this, we measured the influence 
of external tides by considering the rate-of-strain tensor, and espe- 
cially its largest eigenvalue A3, fitted to the flow between 1 and 4 
virial radii around haloes. In detail, we found that for small haloes: 

(i) High formation redshift is related to a low assembly rate at 
late times. Thus, any spatially biased process that suppresses mass 
growth at late times can potentially lead to an age dependent bias 
of halo clustering. 

(ii) Among all measures of environment that we considered, as- 
sembly rate at redshifts z > correlates most strongly with the 
largest eigenvalue A3 of the rate-of-strain tensor, and especially so 
when measured at high redshift, z ~ 1. The correlation of assem- 
bly rate with the density of the environment is weaker. We find 
that haloes with a low assembly rate have suppressed infall or even 
outflow in the region between 1 and 4 virial radii along one di- 
rection. Furthermore, we find that this direction corresponds, with 
little scatter, to the direction to that neighbouring halo that has the 
strongest tidal influence. 

(iii) The strain eigenvalue A3 (probing velocity shear) is 
strongly correlated with the tidal field due to the neighbouring halo 



that exerts the strongest tides. This correlation has been verified us- 
ing a restricted three-body approximation for the tidal field. Thus, 
tidal forces induce a velocity shear in the flow around dark matter 
haloes and can render accretion ineffective by altering the conver- 
gence of the accretion flow onto a halo. 

(iv) The correlation between A3 (probing velocity shear and 
tides) and density grows with time. Thus, haloes in strongly sheared 
regions are much more strongly clustered than average and asso- 
ciated with high density regions onto which they fall. Hence the 
suppression of mass growth is stronger in regions of high density. 

(v) Haloes can undergo mass loss due to the tidal influence of a 
larger halo when their Hill radius becomes smaller than their virial 
radius. This already happens at a distance of 1.5 virial radii of the 
larger halo and will be most severe inside its virial radius. Roughly 
20 per cent of the 2 - 4 x 10 10 h~ 1 M e fulfill the condition for 
mass loss at some point of their history. These haloes are the most 
extreme cases of tidal suppression of assembly. 

(vi) The mismatch between the mass that the Extended Press- 
Schechter formalism would attribute to a halo and the mass that is 
measured for it in the simulation correlates with formation time and 
with the eigenvalue A3 indicating a dependence of this mismatch on 
tides (this is shown in the appendix). 

From the evidence presented above, we conclude that the tidal 
field due to neighbouring massive haloes exerts a shear on the flow 
surrounding small haloes. In the rest-frame of the small halo, this 
is an outward flow along the direction to the more massive halo, 
tidally suppressing the growth of the small halo. This is especially 
effective in filaments since most of the matter is located along the 
direction from which it is least effective for the halo to accrete. 

The proximity to more massive haloes, where the clustering 
is enhanced, leads to a spatial bias of haloes that experience strong 
tides, especially haloes below the non-linear mass, and thus also to 
a bias of the population of "old" and "young" haloes of the same 
mass. We conclude that this is an important factor in the emergence 
of the assembly bias. 
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APPENDIX A: ASSEMBLY BIAS IN THE INITIAL 
DENSITY FIELD 

Al Predicting the assembly bias from the initial density field 

Can we predict the amplitude of the assembly bias directly from 
the linear density field? This requires being able to forecast the 
formation redshift, final mass and local Eulerian density of a halo 
directly from the initial conditions. 



Al.l Halo formation redshift 

What determines the formation redshift of a halo? In Figure |ATI 
we show the Spearman rank correlation coefficient between the for- 
mation redshift zform and the linear overdensity field evaluated at 
the centre of mass of the Lagrangian patch that will form a halo 
at z = 0. The overdensity field is computed after smoothing the 
density field with a top-hat filter of radius R. The rank correlation 
coefficients with Zf or m are plotted against R. Our results show that 
the formation time of a halo of mass M is most strongly determined 
by the value of the linear overdensity smoothed on a scale ~ M/2. 
This is what one trivially expects for spherical perturbations. Tides 
on slightly larger scales (but still internal to the protohalo) also cor- 
relate with Zform but to a lesser degree (not shown in the figure). 
The tidal suppression of halo assembly discussed in the previous 
sections manifests itself in the initial conditions as a residual corre- 
lation between Zt om and 8r for large smoothing scales. The rank 
correlation coefficient is as high as 0.2 even for smoothing radii of 
a few Mpc. The presence of nearby density fluctuations influences 
the formation time of a given halo. 



A1.2 The assembly bias 

If the final halo masses are known, the assembly bias can be per- 
fectly reproduced starting from the initial density field and using 
a simple peak-background split model where small-scale fluctua- 
tions are responsible for halo collapse and large-scale fluctuations 
determine the local density of the halo environment. 

Inspired by Figure IATI we assume that the formation redshift 
Zform of a halo of mass M is fully determined by the "peak height" 
5 m/2- Denoting by D+(z) the linear growth factor of density fluc- 
tuations, we thus solve the condition for collapse of a spherical 
density perturbation 



' <5 M/2 — 5c 



(Al) 



D+{zi) 

where z% is the initial redshift. We find that a value of S c — 2.3 
leads to the best agreement between zf OIU1 and Zf olm obtained from 
the iV-body simulation. 

In a second step, we determine the evolution of the back- 
ground density field 5$ on scales of Rb = 5 h~ Mpc. Small haloes 
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Figure Al. Spearman rank correlation coefficient p s between halo for- 
mation redshift Zf orm and "peak height" of the Lagrangian density field 
smoothed with a real space top hat filter on scale R as a function of the fil- 
ter scale. Results are shown for three mass ranges. Results for the different 
simulations are indicated by different line styles. The shaded regions indi- 
cate the half-mass range M/2 for the corresponding mass bins. This anal- 
ysis reveals: (1) The correlation coefficient peaks as expected for haloes of 
mass M on scales around M/2; and (2) the correlation with the large scale 
density field on several h~ 1 Mpc is significant. 



that survive until z = will preferentially form in regions that are 
underdense on these scales in the initial density field. They will then 
migrate into regions of higher density over time. We determine the 
evolution of the background density using first-order Lagrangian 
perturbation theory. Given the eigenvalues Aj of the initial tidal 
field on scales of Rb,o, the evolution of 85 at the positions of the 
haloes is given before shell-crossing by 

1 



6b(z) = 



UU [!-£>+(*) A,] 



(A2) 



in the Zel'dovich approximation dZel'dovichll 19701) . 

For each halo in the iV-body simulation, we thus obtain an 
environmental density at z — and the redshift of collapse zf or m 
by knowing a priori its final mass and its Lagrangian position. In 
Figure lA2l we show the median formation redshift if orm (in grey) 
as a function of density 65 for three mass bins0 The correspond- 
ing relation is shown also for the quantities directly obtained from 
the TV-body simulation (in black). We find very good agreement 
with the iV-body data suggesting that this simple model suffices to 
predict the correlation between formation redshift and environment 
seen in the simulations if the correct final mass is known. 

Replacing the actual halo masses with those predicted at the 
Lagrangian halo locations by the Extended Press-Schechter (EPS) 
formalism is sufficient to completely erase any dependence of the 
formation redshift on environmental density. The spherical collapse 
model (on which the EPS method is based) tends to overestimate 
halo masses, particularly for those haloes that are located within 

B We have checked that results are not influenced by the small fraction of 
haloes that already experienced shell-crossing on the 5 ft _1 Mpc scale. 



Figure A2. Halo formation time as a function of environmental density on 
scales 5/i _1 Mpc obtained from the TV-body simulations (grey) and ob- 
tained from the initial conditions data using a combination of two spherical 
collapse models for the evolution of the background and the peak (black), 
the mass of the halo is taken from the A-body simulation. Stripped haloes 
are not considered. Solid lines indicate the median, dashed lines the ±l<x 
spread around the median. 



strongly sheared flows. In Figure IA3I we show how the forma- 
tion redshift Zf or m depends on the ratio between the mass A/eps 
and the mass M, where Zf orm and M are determined from the 
A-body simulations. We find that the mismatch increases from a 
median ratio of about 1.5 at low formation redshifts to about 3 at 
hi gh formation redshifts for the considered masses (cf. Figure 3 
in Wa ng et al. 1 1200% . When we exclude the mass-stripped haloes 
(grey lines) the scatter is significantly reduced. The ratio between 
the EPS-predicted mass and the mass found in the A-body simula- 
tions correlates with the largest eigenvalue A3 of the rate-of-strain 
tensor as strongly as with the halo formation redshift. For haloes 
at z = with masses between 2 and 4 x lO lo /i _1 M0 the Spear- 
man rank correlation coefficient between Meps /M and the shear 
eigenvalue A3 is 0.29. It is 0.27 between Meps/M and Zf orm . This 
provides further evidence that tidally suppressed accretion is re- 
sponsible for the existence of the assembly bias. The proximity to a 
massive halo is reflected in the correlation coefficient with density 
obtained with the top-hat filter at z=0. We find that the correlation 
coefficient with Meps /M is 0.22 for density determined on scales 
of 5 /i _1 Mpc, 0.31 for 2 /i _1 Mpc and 0.41 for 0.5 h _1 Mpc. 



A2 Lagrangian assembly bias 

It is interesting to verify whether the assembly bias is already 
present in the spatial distribution of protohaloes in the linear density 
field. The Lagra ngian bias is given by &l ~ b E ~ 1 if 6b is the Eu- 
lerian bias (e.g. lMo & White! 19961) . For M < M, haloes b E < 1. 

Thus, if bE, old > &E, young it follows that ^ b L ,old > b E , young 

and thus the old haloes should have a lower Lagrangian clustering 
amplitude than the young ones. 

To verify this, we compute the two-point correlation of the 
Lagrangian patches of old and young haloes (selected at z = 0) at 
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Figure A3. The ratio of the mass determined in the initial conditions using 
the EPS model to the mass measured in the iV-body simulation as a func- 
tion of halo formation time. Solid lines indicate the median, dashed lines 
the 16th and 84th percentiles. Shown is the relation for haloes that are not 
undergoing tidal mass-loss (grey) and for all haloes (black) in three mass 
bins. 



the initial redshift of our simulations. Figure lA4l shows the correla- 
tion function of the protohaloes that will form objects with masses 
between 2 and 4 x 10 10 /i _1 Mq at z = 0. The corresponding 
correlations at z — have been presented in Figure Q] At z = 
the oldest haloes cluster more strongly than the youngest ones. For 
the Lagrangian patches, as expected, the opposite is true: the old- 
est (youngest) haloes are significantly less (more) clustered than 
the total. This indicates that, already in the linear density field, 
the environments surrounding the Lagrangian patches of young and 
old haloes have different properties. In particular, we find that both 
young and old haloes at these masses form in initially underdense 
(i.e. with a median 5s < 0) regions. However, the Lagrangian 
patches of young protohaloes reside (on average) in more under- 
dense regions which, in a Gaussian random field, are more strongly 
clustered. Using the mean squared overdensities at the halo centres 
as a proxy for their clustering amplitude (i.e. assuming a linear bias 
model independent of halo age for the halo density), B = (<5f ) 1,/2 , 
we find Boid/S a n = 0.8 and B young / B a \\ — 1.2. On the other 
hand, with time, the protohaloes will collapse and flow towards 
denser regions thus changing their clustering properties. At z — 
we find i3oid/B a n = 1.8 and B y0 ung/Baii = 0.33. The reversal 
of the assembly bias is then a simple result of the evolution of the 
halo environmental density. Older haloes have thus always resided 
in more overdense regions than younger haloes. 




r [h -1 Mpc] 

Figure A4. Two-point correlation functions of the Lagrangian centres of 
mass in the initial density field (z = 78.5) for haloes with masses 2 — 4 X 
10 10 h^M-Q at z = (upper panel) and 2 - 4 X 10 11 h~ 1 M Q (lower 
panel). The solid line shows the correlation function for all haloes, while 
open circles and filled squares correspond to the 10 per cent of these haloes 
with lowest and highest formation redshifts, respectively. 



